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Abstract. Tlie problem of anomalous scaling in magnetohydrodynam- 
ics turbulence is considered within the framework of the kinematic ap- 
proximation, in the presence of a large-scale background magnetic field. 
Field theoretic renormalization group methods are applied to the Ka- 
zantsev-Kraichnan model of a passive vector advected by the Gaussian 
velocity field with zero mean and correlation function oc S{t — t')/k'^^'^ . 
Inertial-range anomalous scaling for the tensor pair correlators is estab- 
lished as a consequence of the existence in the corresponding operator 
product expansions of certain "dangerous" composite operators, whose 
negative critical dimensions determine the anomalous exponents. The 
main technical result is the calculation of the anomalous exponents in 
the order of the e expansion (two- loop approximation). 

Keyvirords: Turbulence, Renormalization Group, Operator Product Ex- 
pansion, Anomalous Scaling, Kraichnan's Rapid-Change Model. 

1 Introduction 

Much attention has been paid recently to a simple model of the passive advection 
of a scalar quantity by a Gaussian short-correlated velocity field, introduced first 
by Obukhov [5] and Kraichnan [7]. The structure functions of the field in this 
model exhibit anomalous scaling behavior, and the corresponding anomalous 
exponents can be calculated within regular expansions in a small parameter. 

Effects of intermittency and anomalous scaling are even more important for 
vector fields. In particular, the large-scale intense anisotropic magnetic fields co- 
exist with small-scale turbulent activity in solar wind, see e.g. jlOj and references 
therein. 

In this communication, we discuss the anomalous scaling of magnetic fields 
in the presence of large-scale anisotropy within the framework of the kinematic 
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Kazantsev-Kraichnan model, using the field theoretic methods of renormaliza- 
tion group and operator product expansion. We extend the one-loop results 
derived in |5j to the two-loop order of the e-expansion. 

1.1 Kinematic MHD Kazantsev-Kraichnan Model 

In the presence of a mean component 0" (actually supposed to be varying on a 
very large scale L, the largest one in our problem) the kinematic MHD equations, 
describing the evolution of the fluctuating part 0—9 (x) of the magnetic field, are 

m 

dtO^ + vd9i = e- dv, + e° ■ dv, + uod'^e^, i = i, . . . , d, (i) 

where the term 6° ■ dvi = fi effectively plays the same role as an external force, 
driving the system, with correlator 

{h{x)f,(x'))=5{t-t')C,,{mr), (2) 

where C is some function finite at r = and decaying for r — )■ oo and m = 1/L 
is the reciprocal of the integral turbulence scale. Here and below x = {t,x.} , d = 
{di = d/dxi} , 9^ = didi = Z\ is the Laplace operator, d is the dimensionality of 
X space, v(a;) is the velocity field. Both v and are divergent-free (solenoidal) 
vector fields: diVi = diOi = Q. 

In the real problem, v obeys the NS equation with the additional Lorentz 
force term oc {d x 6) x 0, which describes the effects of the magnetic field on 
the velocity field. The framework of our analysis is the kinematic MHD problem, 
where the reaction of the magnetic field on the velocity field is neglected. We 
assume that at the initial stages 9 is weak and does not affect the motions of the 
conducting fluid: it becomes then a natural assumption to consider the dynamics 
linear in the magnetic field strength. 

More precisely, we shall consider a simplified model, in which v(a;) is a Gaus- 
sian random field, homogeneous, isotropic and (5-correlated in time, with zero 
mean and covariance 

{v.{x)v,ix')) = D,^-^ J dkP,,(k)fc-'^- . e[^'-(---')], (3) 

where Pij(i<.) = 5ij — kikj/k"^ is the transverse projector, k is the momentum, 
k = |k| , Dq is an amplitude factor, d is dimensionality of the x space and e is 
a free parameter. The IR regularization is provided by the cutoff in the integral 
from below at fc = to cx 1/L. The case of anisotropic velocity ensemble was 
studied in [6]. 

2 Field Theoretic Formulation 

This stochastic problem is equivalent to the field theoretic model of the set of 
three fields = {0, 0', with action functional [TT| : 

Si^) = 0'Dg0'/2 + e'^ [-dt9, + vqAO, - dkM - v^Ok)] - (4) 
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where the first four terms represent the De Dominicis-Janssen-type action for the 
stochastic problem UHl]) at fixed v, and the last term represents the Gaussian 
averaging over v. Dq = (ff) and Dy — (vv) are the correlators ^ and ([3]) 
respectively, the required integrations over x = {t, x} and summations over the 
vector indices are understood. 

The diagrams of the perturbation theory are constructed of the four elements. 
In the a;,k representation the factor i[kaSbc — Pb^ac] corresponds to the vertex, 
and the lines vv, 99 and 99' correspond to the bare propagators 

_ P.,(k) _ P.,-(k) ^ C.,(k) 

-J j- fc<i + '! ' fr- -iuj + uk^^r j- aj2_|.j^2fc4, (^5j 

where Cij(k) is the Fourier transform of the function from Eq. 

The UV divergences manifest themselves as poles in e in the diagrams. For 
the complete elimination of these divergences it is sufficient to perform the multi- 
plicative renormalization of the parameters vq and go with the only independent 
renormalization constant (see Ref. [I]): 

i^Q = vZy, go = gii^Zg, Z'^ = Z^. (6) 

The exact response function Gij = {9i9'j) satisfies the standard Dyson equa- 
tion with just one self-energy diagram, and therefore one can obtain an exact 
expression for the renormalization constant Z^, in the MS scheme: 

^. = l-^i^-^, u = gSal{2^)^ (7) 
where Sd is the area of the unit sphere in d-dimensional space. 

2.1 RG Equations for Composite Operators 

The basic RG equation for a multiplicatively renormalizable quantity F = Zp ■ 
Fft (correlation function, composite operator etc) has the form 

[Org + 1f]Fr = 0, Drg = D^+ pdu - juD,, = xd/dx, (8) 

where RG functions /3 and 7 (anomalous dimension) are defined as 

/3 = = In Zf VZp , (9) 

where Z?^ is the operation at fixed bare parameters. 

From the analysis of RG functions it follows, that the RG equations ([5]) 
possess an IR stable positive fixed point : 

= ^e, /3(uO = 0, P'{u,) > 0. (10) 

The value of anomalous dimension 71/ (w) at fixed point is 



7* = 7^(u*) = e. 



(11) 
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This fact implies that correlation functions of this model exhibit scaling be- 
havior; the corresponding critical dimensions A[F] = Ap can be calculated as 
series in e. For the basic fields and quantities the dimensions are found exactly 

Ag = -l + e/2, Ag,^d+l- e/2, A^ = I - e (12) 

(there is no corrections of order and higher, this is a consequence of the exact 
equality 71/ (u*) — e). 

Let G{r) = {Fi{x)F2{x')) be a single-time two-point quantity; for example, 
the pair correlation function of the primary fields <l> — |0, 0', or some multi- 
plicatively renormalizable composite operators. The solution of the RG equation 
gives: 

G{r)'^vf'A'^'='{Ar)-'^''i{mr), (13) 

where the canonical dimensions d^, da and the critical dimension Aq of the 
function G{r) are equal to the sums of the corresponding dimensions of the 
quantities Fi. 

This representation describes the behavior of the correlation functions for 
ylr ^ 1 and any fixed value of mr. The inertial range = I <^ r <^ L = 
corresponds to the additional condition mr <^ 1, the form of the functions ^{mr) 
in the interval mr <C 1 is studied using the operator product expansion (OPE). 

3 Operator Product Expansion 

According to the OPE, the single-time product Fi{x)F2{x') of two renormalized 
operators has the form 

Fi{x)F2{x') =J2Cc.{r)F^{^,t), (14) 

a 

where x = (x + x')/2 = const, r = x — x' — > 0, the functions Ca are the Wilson 
coefficients regular in m^ and Fa are all possible renormalized local composite 
operators allowed by symmetry, with definite critical dimensions Z\q,. 

The renormalized correlator {Fi{x)F2{x')) is obtained by averaging Eq. (|14p 
with the weight expS'j^; hence the desired asymptotics for the correlator 
{Fi{x)F2{x')) is the sum, in which the operator possessing the minimal dimen- 
sion gives the leading term: 

^[mr) = const ■ (mr)'^™'". (15) 

The feature typical to the models describing turbulence is the existence of 
composite operators with negative critical dimensions, for example, critical di- 
mension of the field 9i is Ag = (—1-1- e/2); their contributions in the OPE lead to 
singular behavior of the scaling functions at mr — ;> 0, that is, to the anomalous 
scaling. The operators with minimal Ap are those involving the maximal possi- 
ble number of fields 9 and the minimal possible number of derivatives. Therefore 
the needed operators are tensors, constructed from the fields 6i themselves: 



F^i = e,,---e,,{eAY,n = i + 2p. 



(16) 
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The critical dimension of any multiplicatively renormalizable quantity F — Zp ■ 
Fn \s Ap = df'p + A^df, + 7^,. Then for the operator F„/ we obtain 

Ap^,=n{-l + e/2) + -i*p^^^, (17) 

and the leading asymptotic term of the correlator {Fni{x)Fpq{x')) in the jth 
anisotropic sector has the form 

{Fniix)Fpq{x')) cx (Ar)' "^^r^i - '^'^p, (rnr)'^''"+p.^ . (18) 

Thus one has to calculate critical dimensions Ap^^ of the operators F^i- 

4 Scalarization of the Diagrams 

The operator Fni — 9i-^ ■ ■ ■ {OiOiY, n — I + 2p is renormalized multiplicatively, 
Fni = Zni ■ F^i, and the renormalization constants Z„; — Zni{g,e,d) are deter- 
mined by the requirement that the 1-irreducible correlation function 

(F„^(x)0(xi) • • • 0(x„))^__ = Z-^^ {F,u{x)e{x^) ■ ■ ■ 0(x„))i__ 

= Z^i^rniix;xi,...,Xn) (19) 

be UV finite in renormalized theory, i.e., have no poles in e when expressed in 
renormalized variables ©• 

Below we present, along with respective symmetry coefficients, all the dia- 
grams needed for the two-loop calculation of the function Fni , except for those 
with the self-energy insertions in the 96' lines. 




(^i) {D2) {D^) {D^) 

The contribution of a specific diagram into the functional Fni has the form 

r = ^p,A,;, (20) 

i 

where Pi are known combinatorial coefficients and Ai are certain scalar quantities 
5 Calculation of Anomalous Exponents 

In this section we present the two-loop calculation of the critical dimensions ^p 
of the composite operators Fni, which determine the anomalous exponents in 
expression (jl8|) . We need to extract from the diagrams only the singular parts 
that contain the first-order poles in e. 
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Calculation of Diagram 2^2- The diagram D2 is represented by the integral 

n+oo^,„ /.+oo^,„, ^00 ^oc Pbi(k + q) 



~ J-^ 277 27r L i'^irr Jm 



z(w+w') + (k + q)^ + + A;<*+« 

• (i(A; + g)oi5ic - i{k + q)cSio) ■ {i(k + q)p5je - i{k + q)eSpj) ■ 

■ {ikqSki - ikiSkq) ■ {ikrdmf - ikfSrm) ■ (21) 
The use of transversality of the vertex greatly simplifies our calculations: 

P1V123 = PliP2Sl3 -P3^12) = 0, ^ PmniP + k) ^ Snm- (22) 

Therefore after contraction of the vector indices with standard symmetric 
structures, constructed from 5-symbols, and differentiation over m, which allows 
to single out the first-order pole 1/e explicitly, one obtains: 



de 



(a;2 + 2a;cos^ + l) 
a;2 + l 



1 



- 1 



sm 



sin^6l 



1 

150' 



(23) 



= — , (24) 
36' ^ ' 



_{x'^ + 2xcos6 + l) 

whcac; ,7: = q/rn. All such integrals were calculated analytically for the most 
important physical case d = 2>. 



Calculation of Diagram JD3. In the similar manner, for the diagram one 
obtains: 



= ''Jo 



de 



.V- + 1 



(a;2 + xcos^ + 1) 



{'e if 8^3 656' 



X 48 



5"^"' 75 



dx / d0 



24 



poo pTT 



,^ 4cos^sin^^ 
de— 7r—r+ 



x^ +x cos e + 1 



(25) 
(26) 



x^ + x COS e + 1 



cos^ e sin^ ( 



J_ / V3_ , 24 
'24 ' 
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Calculation of Diagram iDi . For the simplest diagram Di one obtains 

^i=0,A2 = -f (27) 



Diagram D4. The factorized four-ray diagram D4 contains only a second-order 
pole in e and therefore is not needed for the calculation of jp. 
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Anomalous Dimension 7^ ^ . The value of anomalous dimension jp is 

7> = Y.'^PiAiU, + 2p,Aiul), (28) 

i 

where the quantities with bars and with tildes correspond to the one-loop and 
two-loop contributions, respectively. 

Finally, combining ^ with ([26]) and for the anomalous 

dimension of the operator Fni with arbitrary n and I one obtains: 

If/ , n,n , , 2 f2n{n-2) 221(1 + 1) , 



7l^„, =-{^K'^ + 3)-2?(/ + l)].e + 62. 



125 375 



n(n + 3) 3 f 82 V , 

-^3^ + 35 + 75 j - 2) [^"(^ - ^) + ^^(^ + - 



6 Comparison with the Exact Solution 

The exact solution for the pair correlator of the problem ([Ij was derived in [5]; 
see also [5] for a more detailed discussion. In particular, the exponents Co and C21 
describing the scaling behavior in the isotropic and leading anisotropic sectors, 
were derived exactly for any d. Expanding those expressions to the second order 
in e and setting d — 3 gives 

In the RG approach, these exponents should be identified with the anomalous 
dimensions of the operators 9i0i and 9i0j, that is, with "fp^^ and Jp^^ - It is easily 
checked that our expression ([^ is in agreement with (pOjl . 



7 Conclusion 

We have applied the RG and OPE methods to the simple Kazantsev-Kraichnan 
model, which describes the advection of a passive vector by the Gaussian velocity 
field, decorrelated in time and self-similar in space. 

We have shown that the correlation functions of the vector field in the con- 
vective range exhibit anomalous scaling behavior, what is closely related with 
existence in this model of composite operators with negative dimensions. The 
corresponding anomalous exponents have been calculated to the second order of 
the e-expansion (the two- loop approximation). 

It is worth noting that the hierarchy relations between the anisotropic ex- 
ponents [S] persist in the two-loop contributions. It is also worth noting that, 
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in contrast to the scalar case, the two-loop contributions for scalar operators 
have the same sign as the first-order ones, see e.g. (|5n)) for 7^,^^^ = Co- Thus 
the anomalous scaling and the anisotropic hierarchy become even more strongly 
pronounced due to the higher-order contributions of the e-expansion. 

The agreement between the exact exponents for the pair correlation function 
is also established. This fact strongly supports the applicability of the RG tech- 
nique and the e-expansion to the problem of anomalous scaling for finite values 
of e, at least for low-order correlation functions. 
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